Hall constant of strongly correlated electrons on a ladder 
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The Hall constant Rh in a tight-binding model of correlated electrons on a ladder at T — is 
expressed in terms of derivatives of the ground state energy with respect to external magnetic and 
electric fields. This novel method is used for the analysis of the t-J model on finite size ladders. It 
is found that for a single hole Rh is hole-like and close to the semiclassical value, while for two holes 
it can vary with ladder geometry. In odd-leg ladders, Rh behaves quite regularly changing sign as 
a function of doping, the variation being quantitatively close to experimental results in cuprates. 



PACS numbers: 71.27.+a, 72.15.-v, 71.10.Fd 

The Hall response in materials with strongly corre- 
lated electrons remains one of the properties least un- 
derstood theoretically. The subject has been stimulated 
by experiments in superconducting cuprates [Q, reveal- 
ing anomalous doping and temperature dependence of 
the Hall constant Rh{T) in the normal metallic state. 
For instance, it is well established that the Hall effect 
is hole-like, Rh > 0, in materials with a low density of 
holes rih, introduced by doping the reference antiferro- 
magnetic (AFM) insulator. The clearest realization is 
La,2- x Sr x CuOi (LSCO), where the doping x can be di- 
rectly related to the concentration of mobile holes per 
unit cell rih = x and the semiclassical result Rh — 1/ n^eo 
seems to be obeyed at lowest T > T c and at low doping 

Theoretical attempts to calculate the Hall effect in 
models of strongly correlated electrons resulted in quite 
controversial conclusions. Even for weak correlations 
H or for the problem of a single carrier in a Mott- 
Hubbard insulator Q , the analysis of the Hall response 
is fairly involved. In more recent investigations rele- 
vant to cuprates, the dynamical Hall constant Rh{^) 
has been studied within linear response theory for the t-J 
and Hubbard model, analytically by high-w, T expansion 
|^| and numerically via exact-diagonalization studies of 
small systems HQ . The obtained results are quite consis- 
tent for the high-frequency quantity R* H = Rh(^ — * oo), 
showing at high T a transition from a hole-like, R* H > 0, 
to an electron- like, R* H < 0, at a finite crossover n* h ~ 1/3 
||. For the most interesting d.c. limit Rh — Rh(u = 0) 
the majority of results obtained for 2D systems at low 
doping and T — > indicate Rh < (?]] , instead of the ex- 
pected hole-like behavior ||. On the other hand, one of 
the present authors || recently showed that for a single 
hole doped into a 2D AFM at T = the result should be 
the semiclassical one with Rh > 0. 

/From another perspective and stimulated by synthesis 
and experiments on novel cuprates, models of interacting 
electrons on ladder systems have also been extensively 



studied in recent years [|10| . The idea is that ladders with 
a variable number of legs can offer a broader insight into 
the behavior of correlated electrons and thus can lead to 
an understanding of the more challenging 2D systems. 
Again, results for the Hall response Rh(oj) at low doping 
obtained through linear response theory reveal a quite 
puzzling, electron-like, Rh < Jll[] behavior. 

Our aim in this work is to formulate and calculate the 
Hall constant as a ground state (T = 0) property for 
a tight-binding model with a ladder geometry. This is 
possible due to the finite transverse width of the system 
that, in contrast to an infinite 2D (or higher D) system, 
does not require a relaxation mechanism to describe a 
proper transport regime. Such a formulation allows for a 
more transparent calculation of Rh and in particular the 
determination of its sign Q . In the following, we apply 
this method to the t-J model on a ladder. Via a numer- 
ical analysis of small systems, we investigate Rh for few 
holes Nh = 1,2 and a finite concentration of holes rih > 
introduced into an AFM correlated spin background. 

Let us consider the simplest single-band tight-binding 
model of interacting fermions on a ladder geometry with 
M legs in the y direction, L rungs in the x direction. 
Periodic boundary conditions (p.b.c.) are assumed in the 
x direction and unit-cell length ao — 1. To analyze the 
Hall response, the following additional ingredients need 
to be incorporated into the model: 

a) a finite transverse electric field £ y = £ ^ has to be 
taken into account, 

b) a homogeneous magnetic field B perpendicular to the 
ladder, introduced via the Peierls substitution. In the 
Landau gauge A = B(—y, 0), only the phases of hopping 
integrals in the x direction, H x , are modified by a phase 
ip = cqB (% = 1). We also assume that the interaction 
term Hi nt is not influenced neither by A nor by £ . 

c) a steady electric current density j — j x is induced 
in the ground state by piercing a closed ladder in the y 
direction with a flux <E>, modifying the hopping term by 
a phase 9 = eo$>/L. 
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The tight-binding model can be thus written as, 



H — H x + H v + H A + H{, 



M 



,i[(m-m)ip-6] ( t 



( C m,i+l,s Cmi;s + H.c), 



m — 1 is 
M-l 



H v = E E( c ™+i,- c ™- + H - c -^ w 

m—l is 

Ha = A^(m - fh)n mi , 

im 

where A = eoE, fh = (M + l)/2 and Hi nt will be chosen 
later on. 

The idea is to study the ground state energy E(9, A, <p) 
of the system as a function of A, 9, ip in order to evaluate 
the Hall constant Rh, 
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We use the fact that the electric current density j and po- 
larization density P = P y in the ground state |0) can be 
evaluated via derivatives of the energy E(9, A, <p) using 
the Feynman-Hellmann relations, 



e dE 
N~d9' 



p= _eodE 
NdA' 
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where N — LM denotes the number of lattice points. 

In the absence of the magnetic held, (p = 0, we use as 
starting point an equilibrium and nonpolar ground state 
with j = and P = 0. Such a ground state might corre- 
spond to finite values 9$ and A . In ladders one expects 
(for a nondegenerate ground state) A = by symmetry. 
On the other hand, in finite systems, in general we find 
9q ^ 0. This can be considered as a finite size effect, since 
in the thermodynamic limit L — > oo no macroscopic cur- 
rent is expected in the ground state and so 9o — > (for 
a particular study of finite-size scaling of 9q see e.g. Ref. 
p2| ) . Taking a proper starting 9q in the following calcu- 
lations is however crucial for obtaining a sensible result. 

Next, to simulate the Hall effect we analyze systems 
with small but finite current j : ^ imposed by a finite 9, 
magnetic field imposed by a ip ^ and choosing a A/0 
so that the system remains nonpolar, P = 0. Hence we 
have to study the variation of E(9q + 9, A, ip) for small 
9, A, ip. It is enough to consider a Taylor expansion up 
to 3 rd order, simplified by invoking: (i) the symmetry 
of the current operator j = f ^0,(^ = 0) under 
reflection (A — > —A) and, (ii) the reflection antisymme- 
try of the diamagnetic current jft — j± ^{9 ^ 0, ip = 0) 



and the polarization operator P = — 
This leads to: 



0,^ = 0). 



E° e9 P + -E° AA A 2 + E° Av Ap 



E° dA jAp + l -E° eAA 9A 2 + ±E° m e 3 + ■■■. (4) 



The superscript zero indicates derivatives at equilibrium 
(at tp = 0). From Eqs.^Jij) and to leading order, j is 
given by 



3 = %E° 0e ~6, 
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while the Hall field A is set by the condition that P = 
even in the presence of finite 9, ip, 

E° AA A + E Av p + E° eA J<p + E° 0AA 9A = 0. (6) 
Retaining terms linear in ip and 9 we can express A as 



A =A„ 
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We are interested in the second term, i.e. in Aj induced 
by finite 9 and related j. We note also that Eg A can 
be expressed simply as the derivative taken at the origin 



Ao = shifted to A„. Inserting Aj from Eq.(^) and 
from Eq.(|^) into Eq.(g), we obtain the expression 
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This is a central result in this work. The main advan- 
tages of the new approach are: (i) Eq.(||) requires the 
knowledge of only the ground state energy, (ii) the con- 
dition for j = 0, P = in the reference ground state of 
finite size systems is much more transparent. 

Now we will show how this formulation is related to the 
linear response theory for the particular case of ladders 
at T = 0, where Rh is evaluated via the dynamic (in 
general complex) Rh{u) ||, 



R h (uj) 
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<7 a /3 denotes the conductivity tensor evaluated at B ^ 0, 
while a^p at B = 0. {r a p) are stress (kinetic energy) 
tensor components, in ladders nonvanishing only in the 
direction of p.b.c, i.e. for a = /3 = x. 

We are interested in the limit ui — » 0. At T — and 
with p.b.c. in the x direction, a xx describes a dissipa- 
tionless singular conductivity given by 
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where D xx denotes the charge stiffness ||13j,|12j. On the 
other hand, in the y direction the polarizability Xyyi^ ~* 
0) is finite due to open boundaries. Using the relation 
j y = dP/dt, we get from Eq.(pl) 
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£T° (CJ -> 0) = i^Xj, 
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The off-diagonal a yx (w = 0) can be written as 



a yx (0) = -iN dt([P(t)J}) = 
Jo 
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where the ground state |0), excited states |m) as well as 
E$a refer to <p ^ but A = 0. Taking into account that 
Eg a oc ip and inserting relations (10T2|) into Eq.(^) we 
recover the expression (||), provided that E A = 0. The 
equivalence for the case E Atp 7^ can also be obtained if 
one calculates the linear response <r yx in Eq.([l2]) not at 
A = 0, but rather at the proper A = A^. 

/From linear response theory we observe that the ex- 
istence of the simple expression Eq.(|J) is subject to the 
presence of the restricted geometry which implies a finite 
a yx (u> = 0) as well as a finite a xx (uj)a yy (uj) for u> — > 0. 
Both quantities would diverge for T = at the 2D limit 
(M — ► 00), although Rh{u — > 0) is expected to remain 
well defined and bounded ||]. 

Let us first test the method for noninteracting electrons 
on a two-leg (M — 2) ladder. It is here easy to find 
the single-electron eigenenergies e±(k, 0, A, cp), referring 
to the upper and lower band. At T = 0, states in both 
bands are occupied for e ± < cf, and the result follows 
from Eq.§), 
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where n^r are electron densities in both bands and = 
4iS|fc|<fc± cos ^- Note that Eq. ( fl3| ) reduces to plausible 
expressions: a) the semiclassical result Rh = —l/n e eo 
for an empty upper band, n+ = 0, and b) Rh = 
for a filled lower band n~ = 1, where = = 2 — n e 

is the density of holes in the upper band. 

Now, we illustrate the method and expected as well 
as anomalous features of the Hall response in correlated 
systems on a study of the isotropic t-J model (t' = t). 
The interaction term in Hamiltonian ([!]) describes AFM 
exchange interaction between fermionic spins on neigh- 
boring sites j — (to, £), 



Hi, 
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and fermionic operators in the kinetic energy term are re- 
placed by projected ones, forbidding a double occupancy 
of sites. Note that the projection does not influence the 
general formalism, Eqs.(Qj8|). 

To obtain E(9, A, ip) and consequently Rh in finite 
size ladders we employ the Lanczos diagonalization tech- 
nique. For a particular system with given M, L and fixed 



number of holes Nh we first find the energy minimum at 
the equilibrium 8q and then calculate numerically at this 
point derivatives E° A A , E° gg , E° , E° A(p , E° gAA . Then R H 
is evaluated using relations (0,||). 

In Fig. la,b we show results for the dimensionless r# = 
gqRh/N in the case of a single hole Nh = 1 on M = 2 and 
M = 3,4 ladders, respectively, of varying length and as a 
function of J/t. Note that the semiclassical result in this 
case would be = 1. We should note that in general 
we find here 9q ^ 0, 7r. Moreover the influence of A v ^ 
is essential since E gAip differs from E gAip significantly. 
E.g., for M = 2, in the most relevant regime J < t both 
quantities can even be of a different sign. We notice that 
results for different L's are quite consistent. J = is 
a special case with a ferromagnetically polarized ground 
state, S tot = (N-l)/2 and hence r H = 1. Also, as Fig. 1 
shows, for J > Rh is hole- like with r# ^ 1. This means 
that the behavior is very close to the semiclassical one 
||, but the deviation from the latter is finite (although 
smaller for larger M) and seems to persist also for L — » 
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FIG. 1. Dimensionless Hall constant th = eoRu/N vs. J/t 
for a single hole on ladders of different lengths L with: a) two 
legs and, b) three and four legs. 

Results for two holes are presented in Fig. 2. At very 
low doping, rih -C 1, one might expect that in the ther- 
modynamic limit holes behave as independent particles 
so that rjj ~ 1/A/j. This is definitely not the case for 
M = 2, where in the majority of the parameter regime 
we even find Th < 0, consistent with Ref. [jll]]. It seems 
that this phenomenon is related to the existence of a spin 
gap in M — 2 ladders and quite pronounced binding of 
holes into pairs [ ^0[ . Results appear more regular for 
M > 2. Conclusions quite consistent with the semiclas- 
sical rjj ~ 1/2 are obtained for the M — 3 and M = 4 
ladders. It is well known that odd-leg ladders do not 
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show a spin gap JlOlj , so this can serve as an explana- 
tion for the essential difference between the M = 3 and 
M = 2 case. For M = 4, a small spin gap is expected 
in undoped ladder however its effect on rn is not 
visible, at least not for reachable L. 




FIG. 2. r_fj = cqRh/N vs. J/t for two holes on ladders of 
different size L x M. 

In systems with more holes, Nh > 2, and available 
N < 20 we are dealing already with a substantial doping 
In Fig. 3 we show results for the doping dependence 
Rninh)- We concentrate on a more regular three- leg lad- 
der, where even-odd effects in Nh are not pronounced. 
Shown are data for two systems, 5x3 and 6x3. Results 
for both systems are in general quite consistent, with de- 
viations appearing only in the regime ~ 0.2 where in 
particular the value for Nh — 4 on a 6 x 3 system appears 
to be irregular, probably a finite-size effect related e.g. to 
a close-shell configuration. 
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FIG. 3. Hall constant eoRn vs. hole doping n/, for the 
three-leg ladder and J = 0.4 t. Results are shown for L = 5 
and L = 6 systems, where the full line serves as a guide to 
the eye. Dashed line represents the R* H result [5]. Shown are 
also experimental results for LSCO taken from Ref. [2]. 

The interpretation of Fig. 3 is straightforward in two 



regimes. For a nearly empty band < 1, n e — 1 — nh ft 
we recover the semiclassical result, e^Rn = — l/n e . 
Analogous, but only approximate, is the hole-like behav- 
ior for low doping nh <C 1 where eoRn ~ l/ n ^- The 
behavior is very asymmetric between the hole and the 
electron side. The change from a hole-like Rh > 
to an electron-like Rh < appears (with the largest 
scattering of results in this regime) at rit ~ 0.27. Out- 
value for the crossover n* h is close to the crossover in 
R* H = l/4n h - 1/(1 - n^) + 3/4 (in 2D and T -> oo) 
at n* h ~ 1/3 H . In spite of similar values for n* h and a 
quantitative agreement for nh > 0.3, R* H (also plotted in 
Fig. 3) deviates at low doping values by factor of 4 from 
the semiclassical result. 

Although we are dealing with a ladder system, we ex- 
pect that results for odd-leg ladders would be very analo- 
gous to 2D systems. It is therefore not surprising that our 
results for Rn(nh) are both qualitatively as well quan- 
titatively close to experimental ones for LSCO (doping 
range < x < 0.35), where values shown in Fig. 3 are 
taken at T — 100 K We note that experimentally 
the crossover appears at x ~ 0.3 and at low doping, data 
are consistent with the semiclassical Rh ~ l/n^eo- 

In conclusion, we have introduced a novel method 
which allows the evaluation of the d.c. Hall constant 
at T = in correlated systems with a ladder geometry 
solely from the ground state energy. Since the behavior 
of ladder systems is in many respects analogous to 2D 
systems the method can be used to approach the anoma- 
lous and theoretically controvertial Rh in cuprates. Our 
numerical results emerging from odd-leg ladders are in- 
deed surprisingly close to experiments in LSCO. 
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